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ABSTRACT 



Context. The most promising mechanism acting towards damping the kink oscillations of coronal loops is resonant absorption. In 
this context most of previous studies neglected the effect of the obvious equilibrium flow along magnetic field lines. The flows are in 
general sub-Alfvenic and hence comparatively slow. 

Aims. Here we investigate the effect of an equilibrium flow on the resonant absorption of linear kink MHD waves in a cylindrical 
magnetic flux tube with the aim of determining the changes in the frequency of the forward and backward propagating waves and in 
the modification of the damping times due to the flow. 

Methods. A loop model with both the density and the longitudinal flow changing in the radial direction is considered. We use the thin 
tube thin boundary (TTTB) approximation in order to calculate the damping rates. The full resistive eigenvalue problem is also solved 
without assuming the TTTB approximation. 

Results. Using the small ratio of flow and Alfven speeds we derive simple analytical expressions to the damping rate. The analytical 
expressions are in good agreement with the resistive eigenmode calculations. 

Conclusions. Under typical coronal conditions the effect of the flow on the damped kink oscillations is small when the characteristic 
scale of the density layer is similar or smaller than the characteristic width of the velocity layer. However, in the opposite situation the 
damping rates can be significantly altered, specially for the backward propagating wave which is undamped while the forward wave 
is overdamped. 
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CN 1. Introduction 



o 



Dynamism is probably one of the most fundamental charac- 
teristics of space and solar plasmas, flows of plasma being 
observed in high resolution on almost all temporal and spa- 
tial scales. Flows are ubiquitous in active region loops and the 
measurements of their velocities have been provided by instru- 
ments like SoHO (see Brekke et al., 1997; Winebarg er et al.L 
120021) . TRACE (see lWinebarger et al . 2001) and morerecently 
Hinode (see for example Chae et al., .2008,; Ofman & Wand . 
120081 iTerradas et all l2008bh . In general the flow speeds are 
small, and in most of the observations they are sub-Alfvenic, 
typically less that 10% of the Alfven speed. The bulk motions 
are observed along magnetic field lines which outline coronal 
loops. These flows could be generated by some catastrophic 
cooling of coronal loops or are related to some siphon mecha- 
nism arising due to the difference in pressure at the loop foot- 
points. Since longitudinal steady flows carrying momentum and 
providing additional inertia are present in coronal loops it is nec- 
essary to study their effects on the transve rse oscillations ob- 
served in these structur e s (see for example lAschwanden et al.l 
Il999t iNakariakov et al.L 11999'). An effect which is of obvious 
importance for magnetohydrodynamic (MHD) wave theory of 
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loops' dynamics is how the period and the damping time are 
modified by the flow. 

Here we are interested in the damping of the fundamen- 
tal kink mode due to resonant absorption, based on the trans- 
fer of energy from a global MHD wave to local resonant 
Alfven waves, and how the efficiency of the mechanism is 
altered by a stationary flow. In the past, the influence of 
a velocity shear on this process due to a longitudinal flow 
has be ei i stud ied by Hollweg et al. ( 1990); Peredo & Tataronis] 
(ll990j) J^derman & Goossens. (1995); Erdelvi et al. (1995^ ; 
Erd elvi & GoossensI (Il996h; iTirrvetal l (Il998l). M ore re- 
cendv. 'Andries et al.' ('20001); lAndries & GoossensI (1200 lb ; 
Erdelvi & Tarovan (2003a b), have investigated in detail reso- 
nant flow instabilities which can occur for velocity shears signifi- 
cantly below the Kelvin-Helmholtz (KH) threshold. These insta- 
bilities are produced when the frequency of the forward prop- 
agating wave (propagating in the direction of the flow) shifts 
into the Doppler shifted continuum of the backward propagat- 
ing wave. Under these conditions the mode becomes unstable 
and the flow acts at the resonant layer as an energy source. In 
most of the aforementioned studies it has been assumed that the 
wavelength is shorter or similar to the tube radius. This is not 
the case for standing kink oscillations in coronal loops which 
are precisely in the opposite regime, i.e., where the thin tube 
(XT) approximation is applicable. 
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In the present paper we extend the previous studies about res- 
onant absorption in the presence of flow to the situation where 
the TT approximation is valid. We start by reviewing the prop- 
erties of kink MHD waves in a homogeneous tube with an axial 
flow and study the nature of the waves, that change from be- 
ing trapped to leaky and eventually become KH-unstable. Then 
we consider a non-uniform tube whose oscillation is damped 
by resonant absorption and investigate how the damping time 
is modified by the flow. In the following analysis three different 
approaches are i mplemented. Firstly, w e use existing theoretical 
work, mainly bv lGoossens et al. l (fT992l) . to calculate the changes 
in the period and damping rates induced by the longitudinal 
flow using the thin tube thin boundary approximation (TTTB). 
Secondly, under the TTTB assumption we derive a linear ana- 
lytical approximation to the damping rate and thirdly, solve the 
full resistive eigenvalue problem without the TTTB assumption. 
Reassuringly, we find that the three methods lead to essentially 
the same results. 



2. Basic features: uniform tube 

We consider what we can call the standard loop model, a cylin- 
drical axi-symmetric flux tube of radius R with a constant ax- 
ial magnetic field Bq and with a density contrast of pi/pe where 
the indices "i" and "e" describe quantities inside and outside the 
loop, respectively. Inside the loop there is an axial flow denoted 
by V,. For simplicity we assume that there is no flow outside the 
tube so that = 0. We start by recalling the analytical results 
obtained for a uniform loop (no transition layer) in the fi - Q 
case. It is well-known that the effect of the flow introduces a 
shift in the frequency of waves and that the known expressions 
for the dispersion relation without flow can be used by simply 
replacing the frequency lo by its Doppler-shifted counterpart, 
Q. - tj - kv, k being the wavenumber along th e tube. The dis- 
persio n relation of MHD waves was derived in Goossens et alj 
(119921) (see also iTerra-Homem e t al.l '20 031: ISoler et all l2008h. 
In the TT approximation {kR 1) Goosse ns et alj ( 19921) found 
that the frequency of the kink MHD wave modified by the flow, 
ojkf, (see their equation Eq. [83]) is given by 



Ukf = k 



PiVj 
pi + Pe 



U)r. 



(1) 



where cjkj is the frequency of the kink wave modified by the 
flow, and fa^cm i s the f requency in the center-of-mass frame as in 
iHollweg etal1([T990l) . 



PiPe 



(Pi + PeV 



.kX 



1/2 



(2) 



In the above equation a>ic is the classic kink frequency of a thin 
static tube and is given by 



C 2 2 \ 1/2 

\ Pi'^A.i+Pe^A ,e 

I Pi + Pe 



with ol>a being the local Alfven frequency, defined as 



(3) 



(4) 



The frequency a>kf of the MHD waves is degenerate with re- 
spect to the azimuthal wave number in the TT approximation. 
Our interest is in the kink waves with azimuthal number m - 1 
since these waves are the ones that move the axis of symmetry 
and the loop as a whole. The plus and minus sign in front of the 



second term of Eq. ([U represent two different waves, one prop- 
agating in the same direction of the flow (forward) and the other 
propagating against the flow (backward). The presence of the 
flow breaks the degeneracy of positive and negative frequencies 
present when v, = v^. = 0. 
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Fig. 1. (a) Real and (b) imaginary part of the frequency of the for- 
ward {oJkf+) and backward (cokf ) propagating waves as a function of 
the internal flow (v,) in units of the internal Alfven speed. The curves 
intersect at the KH critical flow, see dot-dashed line, given by Eq. 
The forward wave becomes leaky before the KHI criteria, since the fre- 
quency of the mode crosses the external Alfven frequency (oi^j), see 
dotted lines, at the value of the flow given by Eq. (|6} (in the TT ap- 
proximation). The dashed lines represent the Doppler shifted Alfven 
frequencies. In this plot L = lOOR, pi/p^ = 3. 



The second term of Eq. ([T]i contains the condition for the 
Kelvin-Helmholtz instability which occurs for any velocity shear 
in absence of a magnetic field. When magnetic fields are present 
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it is straightforward to see that the square root in Eq. ([TJ is neg- 
ative when the flow is faster than the critical value given by 



(5) 



The above condition means that fast flows compared to the in- 
ternal Alfven speed are required for the K elvin- Helrnholtz in- 
stability to occur (see also Chandrasekhail 1 196 it [Ferrari et all 
Il98lh . An equivalent problem has been studied in the con- 
text of propagating tran s verse waves in coronal jets by 
IVasheghani Farahani et alj 12009'). These authors have also 
found that in the observationally determined range of parame- 
ters, the waves do not undergo either to the KHI or the negative 
energy wave instability. 

However, before the Kelvin-Helmholtz instability occurs, it 
may happen that the frequency of the modes is above the external 
cut-ofF frequency ((L>A,e = kvA,e), meaning that the wave becomes 
leaky. The forward propagating wave becomes leaky when the 
following condition is satisfied 



(6) 



Similar to the KH-instability, fast flows are required to gen- 
erate a leaky wave. C ontrary to the static situation studied by 
iGoossens etaP (l2009l) . in the presence of flow an underdense 
loop is not required to have leaky modes when /3 - 0. Solving 
the dispersion relation in the TT approximation, the damping of 
waves due to leakage is derived from the imaginary part of the 
frequency given by 




4 (pi + Pe) CD, 



The decrement is proportional to the square of kR (kg ~ k), 
meaning that leakage is basically unimportant in the TT approx- 
imation since the damping time (td - l/y) is very large. Thus, 
the damping due to MHD radiation in the presence of flow is a 
very inefficient damping mechanism for the transverse loop os- 
cillations, even in the presence of fast flows. 

An example of the dependence of the real (luir) and imagi- 
nary (7) part of the frequency of the modes (fundamental for- 
ward and backward waves) on the flow is shown in Fig. [T]for 
a particular equilibrium configuration (L = 100/?, pi/pe = 3). 
The domains where waves become leaky or when a KHI oc- 
curs are clearly shown. When the frequencies of the two modes 
merge for increasing velocity shear (v,) the system becomes un- 
stable. However, note that the forward wave is always leaky be- 
fore the system becomes KH unstable (compare also Eq. ||6l with 
Eq. Is)). It is important to mention here that we have not consid- 
ered here the Principal Leaky Mode which i s a very pec uliar 
solution of the dispersion relation (see ICallvl Il986l l2003h and 
instead have focused on the modes than are trapped for a the 
static background. 

Once we know the main effects of the flow on the kink MHD 
waves we need to know in which regime of Fig. [T| we can match, 
for example, the observed standing kink oscillations. The ob- 
servations of flows in coronal loops indicate that they are slow, 
therefore hereafter we focus on sub-Alfvenic flows rather than 
the super-Alfvenic flows that might cause leakage and KH in- 
stabilities. We concentrate on the regime v,7v/i, < 0.1, thus ac- 
cording to the previous analysis both the forward and backward 
waves are always trapped. This also prevents the presence of res- 
onant flow instabilities which occur when the frequency of the 



forward propagating wave shifts into the Doppler shifted contin- 
uum of the backward propagating wave. 



3. Waves in a Non-uniform tube 

Now let us consider a tube with a smooth variation of density and 
flow across the loop cross-section. In particular we consider the 
case when p varies from its internal value p, to its external value 
Pe in the interval [R-l/2, R+l/2] and the velocity changes from v,- 
to in the interval [R-l* /2,R + l* /2]. Under such conditions the 
process of resonant absorption takes place and kink oscillations 
in coronal loops will damp efficiently. The reader is referred to 
iGoossensI (|2008) and references therein for a detail review on 
this kink wave damping mechanism. 

As in the previous Section we concentrate on propagating 
waves, the possible excitation of standing waves in the presence 
of flow is discussed later. 



3.1. The TTTB approximation 

In a non-uniform tube the imaginary part of the frequency (of 
the trapped 

propagating modes) is different from zero due to mode con- 
version at the inhomogeneous layer (a> = (Dr + iy). Some time 
ago Goossens et al. ( 1992) derived an expression for the damp- 
ing rate in the thin tube and thin boundary (/ <K R) approxima- 
tion for incompressible MHD waves. For compressible waves in 
a magnetic cylinder, using the loop model considered here, we 
obtain exactly the same expression, given by (see their Eq. [76]) 



(7) r 



2(pi + Pe) OJcm pirA)\^\ fA 

As usual rA is the resonant position that is calculated from 

Q}(rA) = u\{rA), 



(8) 



(9) 



i.e. where there is a match between the Doppler shifted fre- 
quency and the local Alfven frequency. It is assumed that the real 
part of the frequency of the resonantly damped mode is given by 
Eq. ([TJ. From a physical point of view, the condition given by 
Eq. ([9]l means that the eigenmodes resonantly interact with the 
Alfven continuum, which is Doppler shifted as a result of flow. 
The factor A in the denominator of Eq. ([8]l is 



^^-[Q}ir)-ulir)), 



(10) 



which contains a term with the derivative of the flow in the ra- 
dial direction, absent in the static situation, that can increase or 
decrease the value of A. 

Given a particular density and velocity profile, the different 
variables in Eq. ([S) can be evaluated. For simplicity, we use the 
well known sinusoidal profile for the density given by 



p(r) = 



Pi, < r < R - 1/2, 

P± [(1 + I) _ (1 _ I) sin , R-ll2<r<R + Z/2(ll) 

Pe, ' ' r> R + l/2. 



This convenient profile has 

ies about resonant a bsorptio n 

2002; Van Doorssela ere et all [2004 : 
Terradas et aL .2006) . In order to make the mathematical 



been 



used in many stud- 



Ruderman & Roberta 
Arregui et all [200l 
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approach more tractable we also assume that the variation of the 
flow speed is sinusoidal inside the loop layer, i.e.. 



v(r) = f [l 



sin 



0, 



0<r<R-l*/2, 
^^], R-l*/2<r<R + l*/2, 
r>R + l*/2. 



(12) 



The flow is variable over a layer of thickness /* which is not nec- 
essarily equal to the characteristic thickness of the layer I where 
the density is non-constant. In general, in preceding studies it 
has been assumed a discontinuous flow (at the boundary of the 
loop), except by Erdelvi & TarovanI (12003 a '.'b) who considered a 
linear profile to model MHD waves in the Earth's magnetotail. 
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Fig. 2. Forward (f2+), backward (fi-) and local Alfven frequency 
(thick line) as a function of the radial position. The arrows indi- 
cate the location of the resonances. In this plot v,/v/i, = 0.02, 
IjR = 0.1, l*IR = 0.01, L = IQQR and p,/p, = 3. The dotted 
vertical lines mark the limits of the inhomogeneity of the flow. 

For the profiles given by Eqs. (fTTT l and (fT2l i it turns out that 
Eq. (|9]l is a transcendental equation for the resonant position 
va- This equation is solved using standard numerical techniques. 
Depending on the spatial scales of the density and velocity we 
distinguish two different regimes, I* > I and the asymptotic case 
I* «i I. 

The analysis of the fir st situation is rather simple (see also 
iPeredo & Tataronisl Il990l) . since the the forward propagating 
wave has always a single resonant position in the range R < 
rA < R + 1/2, while the resonant position of the backward wave 
is situated in the range R - 1/2 < ta < R- This behaviour is eas- 
ily understood from Fig. |2] where we have plotted the Doppler 
shifted frequencies and the Alfven frequency as a function of the 
radial coordinate. The resonant positions are located at the inter- 
section of with o)^ (see arrows). Note that Fig.|2]also shows 
that if the Alfven frequency is discontinuous (jump in density, 
I - 0) there are no resonances (implying no damping) since 
will never intersect the curve corresponding to w^. 



Once the resonant position is determined \A\r^ is evalu- 
ated and we finally obtain the value of the damping rate y (using 
Eq. [HI). A useful quantity that we can calculate is the the damp- 
ing per period, given by 



P 



\<^kf\ 1 

Irl 2n 



(13) 



In this expression we use the real part of the frequency given 
by Eq. ([1]). In Fig. [3] (see solid lines) t^/P is represented for 
two different values of the characteristic widths of the layers, 
l/R - I* /R = 0.05,0.1, as a function of the internal flow (re- 
call that both the frequency and the damping rate depend on the 
flow). The curves with positive slope correspond to the forward 
waves while the ones with the negative slope represent the back- 
ward wave. Figure [3] also shows that increasing the strength of 
the flow increases the damping per period for the forward wave 
and decreases it f or the backward wave. Thi s is in agreement 
with the results of iPeredo & Tataronisl (Il990h . the shifts of the 
frequency modify the location of the resonant surfaces in such 
a way that one of the natural modes is closer to the resonance 
while the other is further away from the resonance relative to the 
static situation. As a consequence, one mode is damped more ef- 
ficiently than the other Nevertheless, the effect of the flow does 
not significantly change the damping per period for the regime 
considered here (0 < v,7va, < 0.1). Hence the mechanism of 
resonant absorption is very robust in the presence of an internal 
flow. The thicker the layer, the faster the attenuation (c.f., the 
results for l/R - 0.1 with the results for l/R - 0.05), a result 
akeady known for the static case. We also see that the change of 
td/P is quite smooth with respect to v,-. 

In Fig. |4] we have represented the damping per period as a 
function of I* in units of loop radius for two different values of 
l/R. In this plot, /, the characteristic scale of the density tran- 
sition, is fixed (recall we are still in the regime with /* > I). 
For large values of I* compared to / we see that the dependence 
is quite weak with the thickness of the flow profile. The for- 
ward propagating wave has a larger damping per period than the 
backward propagating wave. However, when I* < I the situa- 
tion is reversed. The curves cross and the forward propagating 
wave is attenuated faster than the backward propagating wave, 
indicating that we are at the threshold of a different regime. 

Now let us concentrate on the situation when I* <K I, i.e., 
we investigate the case with a very steep profile for the equi- 
librium flow velocity. This is different to the regime discussed 
earlier (I* > I) in several aspects. In Fig. |5] we have plotted a 
typical example. It is easy to see that the forward wave can have 
now three different resonant positions (see arrows). Apart from 
the resonance inside the inhomogeneous velocity layer, the for- 
ward wave has two additional resonances, one at R - 1/2 < ta < 
R-l*/2 and the second atR + l* /2 < ta < R + l/2. The backward 
wave still has a single resonance situated, as before, inside the 
velocity layer. Although the situation is more complicated, we 
can still understand the role of the resonances. The derivatives of 
Q? with respect to r at the resonance inside the velocity layer, lo- 
cated around r - R (for both the forward and backward waves), 
become very large (in absolute value) for small I*, thus domi- 
nating over the derivative of a>^ (see Eq. fTOl ). This means that 
the factor \A\r^ is large and tending to infinity (for I* I), there- 
fore y will tend to zero, i.e., they will not produce any damping. 
However, the other two resonances of the forward wave still be- 
have as ordinary resonances since the derivative of the flow is 
zero (they are located outside the velocity layer where the flow 
is constant) and in this situation the total damping of the mode 
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Fig. 3. Damping per period as a function of the flow inside the 
loop for the forward (+) and backward (-) propagating waves. 
The solid lines represent the analytical results calculated using 
Eq. ([TJ and Eq. ([8]l. The dashed lines are the approximations of 
the damping per period using Eq. (fTTj i and Eq. ( |28T l. The dots 
represent the full numerical solution of the resistive eigenvalue 
problem. The horizontal dotted lines are the damping per period 
in the static situation. For the curves with IjR - 0.05,0.1 we 
have used I* jR = 0.05,0.1. 



Fig. 5. Forward backward (Q_) and local Alfven frequency 
(see thick line) as a function of the radial position. The arrows 
indicate the location of the resonances. In this plot Vj/v^,- = 0.02, 
IjR = 0.1 and I* jR = 0.0001, L = IQOR and p,/p, = 3. The 
vertical dotted lines mark the limits of the inhomogeneity of the 
flow. 



will be finite due to the combined contribution of the two reso- 
nances. 
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3.2. Linear approximation of frequency and damping time in 
tiie TTTB approximation 

Visual inspection of Fig. [3] shows that the damping per period 
varies smoothly and appears, to very good approximation, to be 
a linear function of vi/vaj. This has motivated us to derive a 
linear approximation of the frequency cok / and the damping rate 
y as function of vi/vaj. Actually, it turns out that 



kvi 



(14) 



is a more convenient variable for obtaining the linear approxi- 
mation. Since 



1 / \l/2 

V2\ Pi I vaj 



(15) 



it follows that vi/vaj <k 1 is equivalent to x <sc 1. 

The first order approximation to (given by Eq. JSJ) is 



0.05 



0.10 

L'/R 



Fig. 4. Damping per period for the forward and backward prop- 
agating waves as a function of the width of the flow profile, /* . 
We use the same notation as in Fig.|3] In this plot Vj/va; - 0.02. 



0.15 0.20 (^cm - (^k, (16) 

so that the Unear approximation to LOk+ (for the forward wave) is 

Pi 



Pi + Pt 



(17) 



In order to derive a linear approximation to y we need some 
intermediate results to be used in Eq. ([8]). The linear approxima- 
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tion to quantity piiQ.j - cjj^ is 



2 2 

2 Pfpe 



iPi +Pe)' 



2 2 



and the linear approximation to p(rA)A can be written as 



(18) 



1 I (pi - Pe) 

4R(pi+p,) 



1 + ^ 



^PiPe 



Pi ~Pe _^ (Pi + Pe) I 



iPi - Pe){Pi + Pe) Pi + Pe iPi ' Pe) I* 



(28) 



p(rA)A = p(rA) 



(i+y) 

In 



(19) 



Here Ip and l^. are the length scales of variation of density p and 
velocity v,, respectively. They are defined as 



1 



1 



p(rA) 
1 

virA) 



dp 
dr 

dv 
dr 



(20) 



(21) 



Note that Ip and /,, are not equal to the width of the non- 
uniform layer of density / nor to the width of the non-uniform 
layer of velocity I* . E.g, for the sinusoidal profile vi/A) - Vijl 
(since we assume that rA - R) and p{rA) - {p, + Pe)/2, the 
resultant characteristic spatial scales are 



_ I Pi + Pe 



I.-'-. 

n 



n pi - Pe 

The quantity y in Eq. ( fT9] l is defined as 

Pi - Pe 



y '■ 



Pi + Pe 



If we assume that /,, ^ L then 



x_ ^ (1 +y) 



lo 



In this case the linear approximation to 1 lipifA) |A|r,) is 



1 



P('-a)|A|,.^ p(rA)(jJ'. 



■{l-y + xJ). 



(22) 



(23) 



(24) 



(25) 



Here 7 is a factor which measures the relative importance of the 
non-uniformity of the flow to that of density and it is defined as 



(26) 



With all terms approximated linearly, the approximation of the 
damping rate, y, becomes 



n Ip ipi - pef 
8p(r^) Pi-pe 



CJk X 



\ + X 



^PiPe 2p; 

iPi - Pe)(Pi + Pe) Pi + Pe 



(27) 



If we repeat the analysis for the backward wave we obtain the 
same expression with a change in the sign in front of x. 

The final expression for the damping rates of the two waves 
(forward and backward propagating) using the sinusoidal profile 
of density and velocity reduces to 



The + sign corresponds to forward and backward propagat- 
ing waves respectively. When the flow is zero, the second 
term inside the curly brackets vanishes and we recover the for- 
mu la of the damping in the static equilibrium (see for exam- 
ple |^d^m^^Roberti,|2002j). In the non-static case since we 
have assumed that Ip = (/ ^ I*) the last two terms inside 
the square brackets are of the same order and we see that the 
imaginary part of the frequency linearly decreases with the flow 
for the forward propagating wave and increases for the backward 
wave. Using the linear approximation to the frequency (Eq. IfTTl ) 
and damping rate (Eq. ||28]| ) it is straight forward to calculate the 
damping per period. The results are represented in Fig. [3] (see 
dashed lines). We see that the analytical approximations agree 
very well with the full solution based on the calculation of 
and the evaluation of Eq. ([T]i and Eq. (O. Note that the depen- 
dence of the damping rate on the spatial variation of the flow 
across the loop boundary, I*, is present in Eq. ( 1281 ). In Fig.|4]we 
have plotted the results using this expression (see dashed lines). 
Again, we find an excellent agreement between the two curves 
for both the forward and backward waves. 

For the regime /* <K / it is possible to derive useful informa- 
tion from the linear approximation. Using the assumption <K Ip 
it is easy to see that the damping rate of the resonance inside the 
velocity layer (for both the forward and backward waves) is pro- 
portional to I*, which means, as we have already anticipated in 
Section 13.11 that the contribution of this resonance to the total 
damping tends to zero (damping time tending to infinity) for 
a purely discontinuous velocity profile. This is the behaviour 
already found in Fig. |4] for the backward wave. Moreover, we 
can estimate the total damping of the two regular resonances of 
the forward wave (see Fig.|5]l by adding the individual damping 
rates. It turns out that the total damping rate for the forward wave 
is simply 



1 I (pi - Pe) 
2R(pi+Pe) 



m X 



1 _ ^ 

(Ok 



hPe 



Pi - Pe 



ipi - Pe){pi + pe) Pi + Pe 



(29) 



This is the asymptotic value for the forward wave when I* tends 
to zero. Note that this is twice the damping rate of the situation 
with a very smooth velocity profile (see Eq. 1281 when I* » I) 
indicating a more efficient attenuation (half the damping time). 

3.3. Beyond the TTTB approximation: fuii resistive 
eigenvalue problem 

The results of the previous Sections are based on the TT approxi- 
mation. It is known that without flows this approximation works 
very well even for thick layers. However it remains to be con- 
firmed whether this assumption is still valid in the presence of 
flows. For this reason, we go beyond the TTTB approximation. 
In this case we solve the full proble m numerically. We follow 
the approach of iTerradas et alj j2006h . To study the quasi-mode 
properties, the eigenvalue problem given by equations (l)-(5) in 
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iTerradas et alj (l2006l) . plus the additional terms due to the flow, 
is solved. A time dependence of the form e'"' is assumed and 
the problem is solved numerically using a code based on finite 
elements. As boundary conditions we impose that the velocity 
components are zero for r — > oo. In practice, the condition is 
applied at r = r^m-, and then it is necessary to check that the 
results do not depend on this parameter On the other hand, at 
r = it is imposed that dvr/dr = 0, i.e., we select the regu- 
lar solution at the origin, while the rest of the variables are ex- 
trapolated. All the variables and the eigenfrequency are assumed 
to be complex numbers, since we are interested in resonantly 
damped modes. We include resistivity to avoid the singular be- 
haviour of the ideal MHD equations at the resonances. The resis- 
tive eigenvalue problem is solved and we obtain the real and the 
imaginary part of the frequency which must be independent of 
the magnetic Reynolds n umber that we use in the computations 
(|Poedts&Kerneiill99ll) . 

The results of the calculations for /* > / are plotted in Fig. [3] 
(shown by dots). The agreement with the analytical calculations, 
using the TTTB approximation, is remarkable. The numerical 
curves almost overlap with the analytical ones. In Fig.|4](shown 
by dots) we represent the damping per period as a function of 
and find the same behaviour as in the analytical expression. 
With these results we are even more confident about the method 
used in Section lTTI and about the analytical expressions derived 
in Section l372l 

For the regime /* <K / the numerical method we are us- 
ing fails since the thinner the layer (in density or velocity) the 
larger the Reynolds number required for the damping time to 
be independent of the dissipation. A method based on the ap- 
plication of the jum p conditions at the res onance or resonances , 
used for example bvlTirrv et all (Il998h or lAndries etaTI (l2000h : 
lAndries & Goossen^ 20011) . is more appropriate but since this is 



not the main focus of this paper it will not be further investigated 
here. 



3.4. The standing wave problem 

The results presented in the previous sections correspond to two 
propagating waves, one propagating in the direction of the flow, 
a)f+, and the other travelling in the opposite direction, In 
general, an initial perturbation will excite these two modes at the 
same time and the system will oscillate with a combination of the 
two frequencies. If the frequencies are real the superposition of 
the two modes (with the same amplitude and phase) will have 
the following form 



COS((L)f+t)+ COS (cjf-t) - 



2 cos|"^-t'"^-'.Uosf"^-:'"^-'r|. (30) 



This equation shows that the loop has an oscillation frequency 

given by (recall that is negative for slow flows) 

modulated by an envelope that oscillates at the beating frequency 



In the absence of flow we have that w/+ 



and 



from Eq. ( |30] | we find the well known pure standing wave solu- 
tion of the static case. For slow flows (the regime suggested by 
the observations), the frequencies of the forward and backward 
waves hardly differ, and the beating frequency is very small (i.e., 
a very large envelope), meaning that to all practical purposes the 
behaviour of the system is equivalent to that of a standing wave. 
When fast flows are considered the motion of the loop is much 
more complex and it is dominated by the beating frequency. 



It must be noted that fact that the amplitude of the oscilla- 
tion is quickly damped with time, as the observations indicate, 
might favour the formation of a standing wave when the flow 
is present. In this situation, when the damping times of the for- 
ward and backward waves are shorter than the beating period 
(lyl » '^/+ ^'^f-^ ) the envelope of the signal is dominated by the 
attenuation due to resonant absorption rather than by the modu- 
lation due to the beating. 



4. Conclusions and Discussion 

We have studied the effect of a longitudinal flow on propagat- 
ing kink oscillations of a coronal loop and their damping, and 
have shown, in agreement with previous studies, that under typ- 
ical coronal conditions a longitudinal flow, which is highly sub- 
Alfvenic, is unable to produce KH-unstable modes. It was also 
demonstrated that leaky modes are generated by fast flows that 
have velocities comparable to the local Alfven velocity. Since 
observations show that flows are at most, 10% of the Alfven 
speed, this means that forward and the backward waves must 
always be trapped in coronal loops. Moreover, the forward wave 
never enters into the Doppler shifted continuum of the backward 
propagating waves (see Fig.[T^) and so there are no resonant flow 
instabilities for slow flows. Although instabilities due to longi- 
tudinal flows are unlikely to occur in coronal loops, other kinds 
of instabilities, for example pro duced by the azimuthal shear of 
the kink mode are possible ( see ITerradas et all l2008at iTerradasl 
l2009UClack & BallaiLl2009h . 

It was demonstrated that the resonant damping mechanism 
due to non-uniform density and flow at the loop boundary is 
not significantly altered by the presence of the flow as long as 
the scale of inhomogeneity of the flow is similar or larger than 
the scale of inhomogeneity of the density. We derived simple 
expressions for the linear approximation to the frequency and 
damping rate as a function of the flow, for forward and backward 
propagating waves in the TTTB limit. These simple formulae 
are very accurate, since they agree very well with the numerical 
calculations of the full resistive eigenvalue problem. The ana- 
lytical expressions will facilitate future seismolog ical applica- 
tions (along the lines of those proposed by lArregui et al.L 120071: 
iGoossens et al l l2008h . since now the damping rate contains the 
velocity flow as an additional parameter. Using these expressions 
we can estimate the differences with respect to the static situa- 
tion. For example, for a loop with flows of v, = O.lv^, and a 
thickness of the layer in density and velocity of 0.05R, the pe- 
riod decreases a 6% and the damping time increases a 14% for 
the forward wave, while for the backward wave the period in- 
creases a 6% and the damping time decreases a 11% compared 
to the purely static equilibrium case. 

A physically peculiar situation takes place when the flow has 
a sharp transition at the loop boundary (in the limit of /* <s: I). 
The backward wave is transformed into an undamped mode even 
in the presence of a non-uniform density transition. Conversely, 
the forward wave is more efficiently damped due to the introduc- 
tion of two new resonances outside the velocity transition layer 

Finally, we must remark that the problem studied in this pa- 
per is an initial value problem where the wavenumber, k, is as- 
sumed to be ral, and we solve for the complex frequency lj. 
Nevertheless, a more convenient description of certain coronal 
loop problems would require to study the boundary value prob- 
lem, where the frequency is prescribed and one solves for the 
complex longitudinal wavenumber. This is will the subject of a 
future work. 
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